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We put forward a method for optimized distillation of 
partly entangled pairs of qubits into a smaller number of more 
entangled pairs by recurrent local unitary operations and pro- 
jections. Optimized distillation is achieved by minimization of 
a cost function with up to 30 real parameters, which is chosen 
to be sensitive to the fidelity and the projection probability 
at each step. We show that in many cases this approach can 
significantly improve the distillation efficiency in comparison 
to the present methods. 

PACS numbers: 03.67.-a, 32.80.Qk, 42.50.Vk, 89.70.+C 



I. INTRODUCTION 

Entanglement distillation or purification denotes 
the extraction of strongly entangled pairs of qubits from 
a larger number of weakly entangled pairs. The objec- 
tive is to share strongly correlated qubits between distant 
parties in order to allow reliable quantum teleportation 
|| or quantum cryptography 0. All methods of entan- 
glement purification require that the two parties perform 
only local operations on their systems and that only clas- 
sical information be exchanged between them, without 
transferring any additional qubits. The possible local 
operations include (Fig. |ï|): (i) unitary transformations 
whereby each party entangles the partides at their dis- 
posal; (ii) non-unitary projections, whereby each party 
measures a portion of their partides, thus projecting the 
rest of the system onto a new state. These projections 
are usually followed by classical communication of the 
measurement results between the parties. Another non- 
unitary operation is filtering @], whereby a pair is with 
a certain probability either discarded or kept after each 
projection, the probability being dependent on the state. 

The principies of the known distillation schemes have 
been surveyed in ^] . Here we focus on recurrence distil- 
lation methods which use two pairs of qubits as input and 
produce a single output pair: (a) The quantum privacy 
amplification (QPA) method [|| requires that the projec- 
tion of the input pairs on any of the four states of the 
Bell basis (|V' ± ) = (0,1, ±1,0)* and 1^) = 

(1,0,0, ±1)*) be larger than 1/2. The QPA method uses 
a sequence of 7r/2 rotations, controllcd-not operations 
and measurements. (b) In the method of Ref. [Q the two- 
qubit input state is supposed to have a projection larger 
than 1/2 on the singlet state This method uses a 



sequence of unilateral tt rotations of the qubits and a bi- 
lateral XOR operation followed by a measurement. The 
output is a pair whose projection on the singlet state is 
larger than that of the input pair. When used recurrently, 
the QPA method [|| often converges faster towards a Bell 
state than the method of U (for the proof of convergence 
of the QPA see §). 
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FIG. 1. Stages of the bàsic entanglement purification 
scheme: (I) The two parties A and B share two partly en- 
tangled pairs Qi and Q2, (II) they apply local unitary trans- 
forms U a and Ub, and, (III) project out one of the pairs, thus 
obtaining the remaining pair in the state g. 

Both of the aforementioned recurrence methods use 
fixed parameters for the state transformations. It is not 
clear whether these methods are optimal for the states 
which are usable as their input, and, if not, how can their 
performance be improved. There are other states for 
which these methods do not work at all. This is demon- 
strated in Fig. 0, which gives the statistics of applica- 
bility of the QPA method for randomly chosen two-qubit 
states. The random two-qubit density matrices were gen- 
erated as in Ref. Q and the statistics was obtained from 
a sample of 10 6 trials ||. For each of the generated den- 
sity matrices it was checked whether it corresponds to an 
entangled (i.e., inseparable) state, and, if so, whether it 
can be purified by the QPA. The results show that about 
37% of all the possible states are entangled and in prin- 
ciple can be used for entanglement distillation. About 



74% of the entangled states have fidelity (see Sec. II B) 
H F > 1/2 and in principle can, after appropriate ma- 
nipulations, be used for the QPA. Only a 12% fraction of 
the inseparable states, which is about 4% of the sample, 
have a diagonal density-matrix element larger than 1/2 
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in the Bell basis, and can therefore be used directly as 
the input of the QPA algorithm. This suggests that there 
is a vast domain in the space of two-qubit states where 
new approaches to the distillation problem can be useful. 
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FIG. 2. Statistics of the entanglement properties of an 
ensemble of randomly chosen two-qubit states. The areas of 
the rectangles (from bottom to top) are proportional to the 
following probabilities: the entire ensemble (All states), the 
probability that a randomly chosen state is entangled, the 
probability that it fulfills the condition F > 1/2, and the 
probability that it is directly usable as an input of the QPA 
scheme (QPA possible). 

Here we deal with entanglement distillation as a prob- 
lem of optimization, aimed at improving the efficiency of 
the recurrence methods and extending the class for which 
they work. In contrast to the previously discussed recur- 
rence methods, we assume that the transformation pa- 
rameters can be chosen at will to our advantage. Further- 
more, as opposed to the ingenious choice of parameters 
in these methods, we do not have to guess their vàlues: 
the optimal choice of parameters for the unitary opera- 
tions and the projections is obtained by minimization of 
an approximate cost function, which represents a tradeoff 
between maximized probability of the conditional mea- 
surement and the best achievable entanglement. These 
principies are similar to the ones previously used for op- 
timized state engineering po[ . However, the present task 
is still nontrivial, since the number of control parameters 
is large (up to 30), and the choice of a cost function is 
far from obvious. 

These issues are analyzed in Sec. ^ where the opti- 
mization procedure is obtained. This procedure is then 
applied in Sec. III to cases where the QPA is inefíicient, 



to input states with fidelity F < 1/2, and to trapped ion 
qubits. The results are discussed in Sec. IV. 



II. OPTIMIZATION PROCEDURE 

Let us assume that we want to prepare, starting from 
two pairs of qubits with density matrices Q\ and g 2 , a 
single pair with a density matrix g whose entanglement 
is larger than that of both q\ and g 2 , using local unitary 
transformations followed by projections (Fig. [!]). Let us 
denote the local unitary transforms of each party by U a 



and Ub, respectively. If the measurements are performed 
on the partides of the second pair, with the detected state 
being \^) 2 , then the first (purified) pair is transformed 
into the state 



Q = p 2^\U a U bei ® Q 2 UIUI\^)2, 

P being the probability of success, 

P = Tr 2 (V\U a U bgi ® g 2 ulul\^) 2 . 
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Provided that each particle is a two-level system, the 
two-particle unitary transformation of each party be- 
longs to the SU(4) group (considering SU(4) instead of 
U(4) means that we omit the unimportant overall phase). 
Such a group has 15 real parameters. The local trans- 
formations on both sides are thus parametrized with 30 
real parameters. Finding the optimum vàlues for these 30 
parameters would enable us to perform the entanglement 
distillation in the most emcient way. For this purpose, we 
must (i) find a proper parametrization of the transforma- 
tions, (ii) choose a function that quantifies the success of 
the distillation, and (iii) have a suitable method for the 
optimization of this function. 



A. Parametrization of the local unitary 
transformations 

To parametrize the SU(4) local unitary transforma- 
tions, we use a modified version of the scheme in Ref. [jïï| , 
in the form of a product of six SU(2) transformations: 

C/=C/ (1 ' 2) (012>12,X12) (3) 

X C/ (2 < 3) (023, 0, X2 3 )t/ (1 ' 3) (013, Vl3, Xis) 

X (034, 0, X34)C/ (2 ' 4) (0 2 4, 0, X24 ) 

xC/(b4)(^ 4)V)l4)Xl4)) 

where the U^'^((f>ij,ijjij,Xij) transforms represent the 
SU(2) rotations between the i-th and j-th states of the 
4-dimcnsional basis, their clements being 
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(4) 



úncpe-™, U^n n) = cos^e-^. 



We consider the basis states to be |<7i,pi), \gi,e 2 ), 
\ei,g 2 ), and |ei,e2), where \gk) and \ek) denote the 
ground and the excited state of the kth system, k = 1,2. 

The local unitary transformations are thus described 
by a 30-dimcnsional vector X, 



X= (012,^2, XÍ2, 023, X23, 

*?3. <3. XÏ3, 0^4, X34, ^24, X24> #4, Mi, X^, 
ib i b b jb b jb i b Ò ±b b 
f>l 2 , Vl2, Xl 2 , 023; X 23 , 013: Vl3> Xl3> ^34; X34> 

024,X24>Ï4^Í4,XÏ 4 ), 



(5) 
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where the indices a and b refer to the two parties (Al- 
ice and Bob, respectively). The distillation protocol of 
Ref m corresponds to the vector 

X! = ~(1, 0,0,0,0,-1, 0,0,0,0,0,0,1,0,0; (6) 
0,0,0,0,0,1,0,0,0,0,0,0,0,0,0) 

(after the randomization yielding the Werner state), 
whereas the QPA protocol || is equivalent to the vec- 
tor 



X 2 
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Thus, both protocols represent specific choices of the 
available transformation parameters. The transforma- 
tions with 15 parameters on each side are the most gen- 
eral possible, but even transformations with fewer de- 
grees of freedom can be suitable for the extremization 
of the distillation efficiency. The number of available 
parameters depends on the particular realization of th e 
qubits and the way we manipulate them (see Sec. III C). 



B. Quantifying the result of the distillation 

The resulting state should be as strongly entangled as 
possible and obtainable with a reasonably high proba- 
bility. The calculation of the probability of success in 
Eq. (H) is easy, but quantifying the entanglement is a 
non-trivial task for which many different measures have 
been suggested. Since finding the extremum of a function 
is time consuming, we prefer a measure of entanglement 
that can be calculated as fast as possible. This has led 
us to choose as our measure the entanglement of forma- 
tion E 0], defined as the least mean entanglement of 
ensembles of pure states realizing the mixed state g (en- 
tanglement of the pure state being the von Neumann en- 
tropy of the reduced one-party density matrix) JT^ . An- 
other convenient measure of the entanglement is fidelity 
(or fully entangled fraction) F, defined as the maximum 
max(e|g|e) taken over all completely entangled states |e) 

IS 

Along with the entanglement, our "cost" function 
should optimize the probability of success. We have ex- 
perimented with the maximization of E and F , and min- 
imization of variously constructed cost functions depend- 
ing on F, E and the success probabilities P. The best 
results have been achieved by means of the cost function 



f c = l/[(F«P + e)(É + e)}, 



(8) 



where F is the largest fidelity of the 4 possible outcomes 
of the measurement, P and E are the corresponding prob- 
ability and entanglement, e is a small parameter (« 10~ 6 ) 
for regularization of the function around P w and E 
w 0, and q is a parameter quantifying the preference for 
larger fidelity or larger probability (typically, q = 12). 
The choice of this function ensures that a large entan- 
glement is achieved with a reasonable probability. By 
manipulating the shape of the cost function (e.g., vary- 
ing q) we can get the resulting state with large entan- 
glement but small probability or vice versa with various 
intermediate possibilities. 

For a comparison of the results of different methods, it 
is useful to estimate the average number of pairs which 
is consumed in order to get one pair with the requircd 
entanglement. Assume first that we have reached the 
goal after nk steps. The joint probability of success in 
all steps is Pfe, which is the product of their individual 
success probabilities. The index k refers to the particu- 
lar sequence of results in the individual steps (different 
sequences of intermediate states can lead to the same 
requircd entanglement). After each step (except the last 
one), the resulting state is taken as the input state for the 
next step. If we have initially Af pairs, then the average 
number of resulting pairs distilled following the sequence 
k would be 



A4 = jt-N, 

2 n k ' 



and the total average number of distilled pairs is 



(9) 



(10) 



where the summation runs over all the sequences k which 
lead to a pair with the required entanglement. The de- 
nominator 2™ fc in Eq. (||) reflects the fact that in each 
step two pairs are consumed to produce one resulting 
pair. From equations (0) and (10) it follows that the to- 
tal number of pairs needed to create one resulting pair 
is, on average 



N 



N_ 
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E 



(11) 



C. Extremization procedure 

We have searched for the extrema of the cost func- 
tions numerically, using the Matlab procedure FMINS. 
This procedure starts from a given point and uses the 
Nelder-Meade simplex search algorithm. Since the ex- 
tremized cost function is generally not convex, the pro- 
cedure finds a local extremum which need not be the 
global one. Thcrcforc, getting a result does not mean 
that we actually found the optimal mcthod for distilla- 
tion. In our computations, we usually begin with several 
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randomly generated starting points and choose the best 
result. In general, we can call it a success if the distilla- 
tion efficiency exceeds that of the methods suggested so 
far |,| . 



III. RESULTS AND APPLICATIONS 

In order to compare our approach to the present meth- 
ods @,|, we first applied the optimization procedure to 
the class of states on which they mostly focus, i.e., the 
Werner states |Q , which are mixtures of the totally en- 
tangled and totally mixed states. In this case the numer- 
ical optimization brought no improvement and the QPA 
method seems to be as efficient as ours for the Werner 
states. By contrast, substantial improvement was found 
for states such that the QPA method either converges rel- 
atively slowly or cannot be used at all || . This refers to 
the states that do not have any diagonal matrix element 
that is larger than 1 /2 in the Bell representation. If onc 
of the diagonal elements in the Bell basis is only slightly 
larger than zero, the QPA convergence may be too slow 
for efficient applications. Let us study these cases in more 
detail by considering characteristic examples. 



50.9% 



0.8% 




0.6 



°- 5 o- 



4 6 
Number of íterations 



FIG. 3. Fidelity as a function of the number of Ítera- 
tions for the optimized scheme (sòlid line) and for the QPA 
(dashed line) with the initial state given by Eq. (^). Joint 
probabilities of obtaining the required states are indicated. 



A. Cases when the QPA is inefficient 

As the first example, consider the state 

Qi = Q2 = §|7)<7l + èhTXVrl + JL±I, (12) 

where I7) = 2 _1 / 2 (0, 1, z, 0) 1 ' and / is the 4x4 identity 
operator. The overlap of this state with the Bell state 
is 17/33 w 0.5152, which is marginally sufficient for 
using the QPA algorithm, and the fidelity defined above 
is F = 0.7518. Let us assume that the aim is to ex- 
ceed F = 0.9, which may be sufficient for application of 



other distillation schemes, e.g., the hashing method [||. 
A direct application of the QPA method would reach this 
value after 9 steps, the joint probability of success in all 
steps being 0.81% (see Fig. |J). This would require an 
average number of 63xl0 3 pairs to get a single output 
pair. On the other hand, our optimization scheme would 
reach the required fidelity in 2 steps with a joint proba- 
bility of 50.87%, so that less than 8 pairs on average are 
consumed to get one output pair, which clearly means a 
substantial improvement of distillation efficiency |15|. 



B. Application to states with fidelity < 1/2 

If the entangled pairs have fidelity F < 1/2, one can- 
not directly use the existing methods |ï],|| to distill the 
entanglement. So far, the only suggested scheme to han- 
dle such pairs would be to transform them first by non- 
unitary operations such as filtering j^] , in order to reach 
fidelity above 1/2. To our knowledge, no explicit for- 
mula for determining the filtering parameters for arbi- 
trary states has been presented. 
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FIG. 4. Purification of inseparable states with F < 1/2. 
The vertical coordinate of the circle of each arrow denotes 
the fidelity of the initial state, whereas that of its tip cor- 
responds to the fidelity of the resulting state after one step 
of optimized distillation. The horizontal coordinate of each 
arrow indicates the probability of success of the distillation 
procedure. The Figure shows 25 examples of the optimized 
distillation procedure for randomly generated initial states. 

Our approach allows for entanglement distillation of 
pairs with F < 1/2 in the same way as for any other en- 
tangled states. To demonstrate this, we have randomly 
generated several density matrices with E > and F < 
1/2, and optimized the local unitary transformations so 
as to reach a state with F > 1/2 (see Fig.|j). We have 
observed that, whereas the success of the optimization 
depends on the value of E in the initial state, it still 
works for all randomly generated states with entangle- 
ment above 5xl0 -4 , allowing inseparable states with F 
< 1/2 to be purified using only unitary transformations 
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and conditional measurements, without filters. Notice 
that starting from a fixed value of F, the distillation can 
exhibit either a large increase of F with a small proba- 
bility or vice versa Jl5|. 

C. Application to trapped ion qubits 

Significant improvement in distillation efHciency is 
achieved by our method not just for rather special, but 
also for genèric, physically important cases. For instance, 
let us consider qubits that are realized by two infernal 
states of trapped ions p6|| . ff two or more ions are 
trapped in a single trap, then the logical functions be- 
tween two qubits are achieved by using an auxiliary in- 
fernal state of each ion and a vibrational mode of the 
collective oscillations. It is assumed that the evolution of 
the iònic states is driven by coherent làser pulses whose 
amplitudes, phases and durations can be controlled. An 
arbitrary unitary transformation of a single qubit can 
be achieved by two resonant làser pulses focused on the 
corresponding ion JÏ6| . The durations (or strengths) of 
the two pulses and the phase of (say) the second pulse 
represent 3 parameters of the transformation. The in- 
teraction between two ions is achieved by a sequence of 
three pulses, whose effect is to flip the sign of the state 
(0,0,0, l) t (i.e., |e 1 ,e 2 }) without changing the other ba- 
sis states. The parameters of these three pulses are fixed 
at properly chosen vàlues so as to ensure that at the end 
of the procedure no auxiliary state remains excited (see 
for details). 
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FIG. 5. Fidelity as a function of the number of iterations 
for the optimized scheme (sòlid lines) and for the QPA scheme 
(dashed lines). Entanglement distillation for two difFerent ini- 
tial states is repeated until the fidelity F = 0.9 is reached. 
The upper tree (starting with F = 0.779) corresponds to a 
decayed singlet (after a decay time jt = 0.25) as the initial 
state. The lower tree (starting with F = 0.691) corresponds 
to the same decayed singlet mixed with a completely mixed 
state as the initial state. Joint probabilities of obtaining the 
required states are indicated for the optimized scheme. 



Let us call this transformation U c , and assume that 
during one step of distillation each party performs only 
one transformation U c , preceded and followed by single- 
qubit rotations. Rotation of all 4 qubits before the U c 
transformation represents 4 x 3 = 12 parameters. Af- 
ter the U c transformation it is sufhcient to rotate only 
the two qubits which are to be measured, as the local 
transformations of the remaining pair do not change the 
entanglement. Thus, we are left with 12 + 2 x 3 = 18 pa- 
rameters to be optimized, their physical meaning being 
the areas and phases of the làser pulses. 

As an example of entanglement decoherence, let us 
take the iònic excitation to be decaying according to the 
(zero-temperature) master equation 

Q = -f [H, g] ~ i (<T+<T-Q + QCT+CT-) +J(T-Q<T+. (13) 

Here the (single-particle) Hamiltonian is H = Tilüq |e)(e|, 
with a + = \e)(g\ and er_ = |g)(e|, where \g) and |e) de- 
note the ground and the excited states, respectively. The 
upper "tree" in Fig. [j| has as its starting point the state 
obtained by the decay of the singlet \ip~) according to 
Eq. (|Ï3|), after a decay time jt = 0.25. This state has 
fidelity F — 0.779. Again, we can assume that the aim is 
to purify the state so as to rcach fidelity of at least 0.9. 
The QPA procedure would achieve this goal in 2 steps, 
consuming on average 7.7 pairs and ending up with a 
state whose fidelity is F = 0.918. Using the optimization 
scheme with 18 parameters, the pairs have a relatively 
large probability to be fully entangled (F = 1.0) after 
the first or second step. We are able to reach full entan- 
glement in this case because the decayed singlet has zero 
probability for both qubits to be in the excited state, as 
opposed to a Werner state. Following the different tra- 
jectòries of our procedure, we find that the mean number 
of pairs consumed before reaching F = 0.9 is only 4.6. 

The lower tree in Fig. || starts from a state obtained 
by mixing the decayed singlet after jt — 0.25 with a fully 
mixed state in a 83% to 17% proportion, as a result of 
additional sources of errors. In this case, even the op- 
timization scheme yields only partially entangled states. 
However, whereas the QPA procedure would need 25.8 
pairs on average to get the required fidelity, the optimiza- 
tion scheme would consume only 15.8 pairs on average for 
this task . 

IV. DISCUSSION 

Our main achievements can be characterizcd as finding 
a straightforward method for efficient distillation of en- 
tanglement, which is particular ly valuable in cases where 
previously suggested methods either do not work or con- 
verge relatively slowly. There is no special requirement 
on the form of the initial states (such as the Werner 
states), except that the states must be entangled. It is 
not even required that the fidelity should be larger than 
1/2. We have seen that this approach allows essential 
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saving of the "raw material" of initial partly entangled 
states. 

Several points must be noted: (i) The optimization 
procedure may end in a local extremum of the cost func- 
tion, which usually requires several trials before the final 
choice of the transformation parametcrs. (ii) By con- 
trast to the previous schemes (ï]||, our method is state- 
dependent: before starting the distillation we have to 
know the initial density matrix. This knowledge is con- 
sistent with the objective of protecting particular cor- 
related two-qubit states, e.g., singlets, from being spoilt. 
The knowledge of the initially spoilt state can be achieved 
either by familiarity with the dissipatio n or e rror dynam- 
ics (e.g., zero-temperature decay - Sec. III C| ), or by state 
reconstruction methods (the problem of density matrix 
reconstruction by local measurements of a pair of two- 
level systems has been studied in detail in Ref. 17 1). Of 
course, in the latter case a portion of the pairs will be 
consumed for the reconstruction measurements, but once 
the density matrix is determined with sufficient precision, 
the distillation scheme can run indefinitely. Note that the 
knowledge of the density matrix would be necessary also 
in the case of the QPA method when F > 1/2 but no di- 
agonal density matrix element in the Bell representation 
is larger than 1/2: the state must be proper ly rotated be- 
fore the QPA method itself is used, and the parameters 
of the rotation would depend on the initial state. 
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